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ABSTRACT 

Most exoplanets detected by direct imaging so far have been characterized by relatively hot (>1000 K) 
and cloudy atmospheres. A surprising feature in some of their atmospheres has been a distinct lack of 
methane, possibly implying non-equilibrium chemistry. Recently, we reported the discovery of a planetary 
companion to the Sun-like star GJ 504 using Subaru/HiClAO within the SEEDS survey. The planet 
is substantially colder (<600 K) than previously imaged planets, and has indications of fewer clouds, 
which implies that it represents a new class of planetary atmospheres with expected similarities to late T- 
type brown dwarfs in the same temperature range. If so, one might also expect the presence of significant 
methane absorption, which is characteristic of such objects. Here, we report the detection of deep methane 
absorption in the atmosphere of GJ 504 b, using the Spectral Differential Imaging mode of HiCIAO to 
distinguish the absorption feature around 1.6 /im. We also report updated JHK photometry based on 
new iG-band data and a re-analysis of the existing data. The results support the notion that GJ 504 b has 
atmospheric properties distinct from other imaged exoplanets, and will become a useful reference object 
for future planets in the same temperature range. 

Subject headings: planetary systems — stars: solar-type — techniques: photometric 
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1. Introduction 

Several companions to stars in or near the plan- 
etary regime have been discovered through direct 
imaging in recent years (e.g. Marois et al. 2010; 
Lagrange et al. 2010; Carson et al. 2013), which has 
provided the opportunity for characterization of their 
atmospheres. The HR 8799 system in particular 
has been the subject of several studies using spec- 
troscopy and multi-band photometry. These obser- 
vations have provided indications of atmospheric 
properties such as non-equilibrium chemistry (e.g. 
Jansonetal. 2010; Hinz et al. 2010; Bowler et al. 
2010) and thick clouds (e.g. Barman et al. 2011; 
Currie et al. 2011; Skemeretal. 2012). The latter 
implies that there may be a difference in cloud re- 
tention between ~1000 K planets and more mas- 
sive brown dwarfs in the same temperature range. 
Brown dwarfs undergo the so-called L/T transition 
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(e.g. Kirkpatrick 2005) at these temperatures, which is 
thought to indicate a transition from cloudy (L-type) 
to clear (T-type) photopheres (e.g. Leggett et al. 2000; 
Chabrier et al. 2000; Allard et al. 2001). The appar- 
ent difference between planets and brown dwarfs in 
this regard can be theoretically justified through a de- 
pendence of the transitional temperature on surface 
gravity (Metchev & Hillenbrand 2006; Marley et al. 
2012 ). 

Recently, we presented the discovery of GJ 504 b 
(Kuzuhara et al. 2013), a planetary companion to a 
nearby GO-type star, made with direct imaging as part 
of the ‘Strategic Exploration of Exoplanets and Disks 
with Subaru’ program (SEEDS; Tamura 2009). With 
an age of >100 myr and a temperature of less than 
600 K, GJ 504 b allows for study of a new class of 
planetary atmospheres, where the clouds have largely 
disappeared from the photosphere, more clearly re- 
vealing the molecular contents of the atmospheres. 
In this paper, we present new data in intermediate- 
band methane filters within //-band and re-analyze the 
broad-band photometry from Kuzuhara et al. (2013) in 
order to acquire a uniformly reduced data set. In the 
following, we first describe the new observations that 
have been acquired, and the reduction procedure for 
the acquired data, as well as the re-reduction of the 
previous broad-band data. We then summarize the re- 
sults, and compare them to a grid of models in order to 
constrain physical parameters, and finally we discuss 
the broader implications of our results. 

2. Observations and Data Reduction 

Here, we discuss the J and //-band observations 
that have already been presented by Kuzuhara et al. 
(2013), new observations in CH 4 S and C ll\ L, and 
a re-observation in /C s -band. Cl 1,\ S and 677] L 

are intermediate-band methane filters ranging from 
1.486 /rm to 1.628 fim for the former and from 
1.643 fim to 1.788 /rm for the latter. These filters are 
designed to capture both sides of the 1.6 /rm methane 
bandhead expected in cool giant planet atmospheres 
under chemical equilibrium. If strong methane absorp- 
tion is present, there should be significantly more flux 
in the CH 4 S band, which is outside of the methane 
absorption feature, than in the CH 4 L band, which 
is inside of it. The CH 4 S and CH 4 L images were 
acquired on 15 May 2012 using the HiCIAO cam- 
era (Suzuki et al. 2010) with the A0188 adaptive op- 
tics system (Hayano et al. 2008) in the 2-channel SDI 
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(simultaneous spectral differential imaging) mode, in 
which the incoming beam is split into two beams us- 
ing a Wollaston prism. One beam is passed through 
the CH 4 S filter and one through the C 1 1,\ L filter, and 
the resulting two images are placed side by side on the 
detector. As with all the other observations of GJ 504, 
the Angular Differential Imaging (ADI; Marois et al. 
2006) mode was used, where the image rotator is used 
to keep the orientation of the telescope pupil fixed on 
the detector (pupil tracking). An exposure time of 20 
seconds was used for each frame, without co-adding. 
During the continuous sequence of 106 minutes (in- 
cluding overheads), 219 frames were acquired with 
these settings, spanning a field rotation of 94.4°. In 
addition, non-saturated images were acquired inter- 
spersed with the ADI sequence (50 frames in total) for 
photometric referencing. These used 1.5 s exposures 
(without co-adding) and a neutral density (ND) filter. 
The ND filter has a transmission of 0.859% across the 
broad //-band, but varies across wavelength such that 
the transmission is 0.714% in CH 4 S and 0.929% in 
CH^L. The /\ s -band observation was acquired on 11 
May 2012 in the same way as for the previous K s - 
band measurement of Kuzuhara et al. (2013) from 28 
Feb 2012, using HiCIAO in the ADI mode. Each co- 
added exposure consisted of 10 co-adds with individ- 
ual integration times of 10.0 s. A total of 38 co-added 
exposures were acquired, spanning a field rotation of 
68.1°. All observations used a circular occulting mask 
with a diameter of 0'/4. 

Since the Wollaston prism splits up the beam based 
on polarization state, there is in principle a risk that ei- 
ther instrumental polarization or polarizated light from 
the planet itself could impact the apparent photometry 
in the SD1 imaging. However, these effects are small 
(note that the planetary flux is emitted light rather than 
reflected light), particularly as the field rotation during 
the observation was more than 90°, which means that 
any polarization effect should be efficiently averaged 
out across the observation. Hence, we treat these ef- 
fects as negligible in this analysis. 

Data reduction was performed using the ACORN S- 
ADI pipeline in a uniform manner, described in 
Brandt et al. (2013). In addition to the CH^S, CH^L, 
and ir s -band data, the J and H - band data from 
Kuzuhara et al. (2013) were re-reduced as well, since 
the discovery paper used a different pipeline and dif- 
ferent reduction parameters and we wanted a uniform 
procedure. We also wished to examine whether un- 
certainties in the centroiding could have led to slightly 


higher photometric uncertainties than had previously 
been accounted for. Consequently, particular attention 
was committed to the registration procedure, which 
was applied after standard de-striping and flat fielding 
of the data. Errors in image centroiding and registra- 
tion propagate into the derived companion photometry, 
and always result in a lower measured companion flux. 
An incorrect absolute centroid will dilute the compan- 
ion flux in a de-rotated and collapsed image. The 
impact of this effect scales with the amount of field 
rotation during the observation. Errors in image regis- 
tration reduce the measured companion photometry as 
shown in Figure 6 of Brandt et al. (2013); for exam- 
ple, root-mean-square errors of one HiCIAO pixel in 
both the horizontal and vertical directions reduce the 
measured near-IR photometry by ~20%. 

The images of GJ 504 were nearly all taken with 
an opaque mask covering the central star, making reg- 
istration and centroiding difficult: a positional shift by 
the star relative to the mask does not simply appear as a 
positional shift in the image on the detector. We there- 
fore started out by not attempting to register the raw 
masked images, nor attempting to fit an absolute cen- 
troid directly. Instead, we initially assumed the frame- 
to-frame pointing stability to be perfect, estimated the 
centroid by eye, and then performed LOCI (Locally 
Optimized Combination of Images; Lafreniere et al. 
2007). We then de-rotated the image sequence around 
many possible absolute centroids and co-added each 
resulting sequence. In this way, it becomes possible 
to find the centroid offset that maximizes GJ 504 b’s 
signal-to-noise ratio in the final, reduced frame. 

Our method partially compensates for a real, sta- 
ble drift of the stellar centroid, by also fitting for the 
drift speed and direction that maximizes the compan- 
ion signal-to-noise ratio. Typical centroid drifts over a 
long observation are generally not more than 1-2 pix- 
els (Brandt et al. 2013), with a ~2-3 pixel drift ob- 
served over a multi-hour observation at very high air- 
mass (Thalmann et al. 2011). We therefore expect the 
uncorrected bias from a stable drift to be much less 
than the 20% that would have been introduced by po- 
sitional fluctuations of 1 pixel, as mentioned above 
(Brandt et al. 2013). Likewise, the actual frame-to- 
frame positional fluctuations in SEEDS images are 
small, <0.3 pixels (Brandt et al. 2013) and will have 
little effect on our photometry. 

Major variations in atmospheric conditions and 
AO performance over the course of an observing run 
present another photometric difficulty. While it is pos- 
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sible to estimate the data quality from the diffraction 
pattern on masked frames, HiCIAO does not introduce 
an unsaturated central peak (as would a semitranspar- 
ent mask) for a more precise calibration. Unsaturated 
photometric frames were sometimes taken only before 
and after an observing sequence and sometimes inter- 
spersed within a sequence. While most observations 
had stable conditions and consistent photometry across 
the sequences, the May 2011 if -band data presented 
a more difficult case. The unsaturated photometric 
frames taken before and after the sequence of masked 
science frames showed a remarkable uniformity in AO 
performance, and the science data themselves seemed 
to show no major AO failures or changes in observing 
conditions. However, using only the first half of the 
sequence yielded a planet flux ~50% higher than that 
derived from the latter half of the sequence. This dis- 
crepancy corresponds to a 2cr significance level. None 
of the other data sets showed any differences of this 
magnitude for companion photometry in subgroups of 
science frames. 

All frames were subjected to a LOCI-based PSF 
(Point Spread Function) modeling and subtraction 
with standard ACORNS-AD1 parameters (AVwn.vr = 
0.6, iVArea = 200). We also tried varying the parame- 
ters, but since the final image is not strongly dominated 
by PSF residuals at the location of the companion, this 
had a negligible effect. All the PSF-subtracted data 
were de-rotated to a common field orientation and 
co-added to form the final image of each filter. The 
primary flux (rescaled from the neutral density filter 
images in the case of CH 4 S and 677,1 L) as well as 
the secondary flux were calculated in circular aper- 
tures with a diameter equal to the 75% of the FWHM. 
Partial subtraction as a function of angular separation 
imposed by LOCI was accounted for when determin- 
ing the contrasts. 

3. Results and Discussion 

The final C / I 4 S and C 1 14 L images resulting from 
the reduction are shown in Fig. 1, and the photomet- 
ric values are shown in Table 1. It can be immedi- 
ately seen from the images that the planet is bright 
in CH 4 S but disappears in CH 4 L. This corresponds 
exactly to the behavior expected from an atmosphere 
with methane absorption, given that CH 4 L encom- 
passes the methane aborption feature, whereas CH 4 S 
represents the continuum flux. Apparent magnitudes 
for the companion were calculated as the measured 


contrast to the primary plus the magnitude of the pri- 
mary in the given filter. Since there are no literature 
values for the primary’s CH 4 S and CH 4 L magni- 
tudes, we translated the //-band magnitude into these 
bands by relating the stellar flux to the Vega flux in 
the relevant passbands, assuming blackbody radiation 
with a T eff of 6234 K for GJ 504 (GOV) and 9602 K for 
Vega (A0V). As expected, the corrections were very 
small, with C H 4 S— H = 0.02 mag and CH 4 L — H = 
—0.01 mag. Apparent magnitudes were translated 
into absolute magnitudes using the parallax distance of 
17.56±0.08 pc from Hipparcos (van Leeuwen 2007). 

The updated J and H photometry is ~0.1 mag 
brighter than the previous photometry. This is consis- 
tent within the error bars, but we adopt the new pho- 
tometry as it more accurately represents the full set of 
flux losses possible in the reduction procedure. Our 
Ft's -band photometry based on the newer (May 2012) 
data point alone is ~0.3 mag fainter than previously, 
which is again consistent given the larger scatter in 
the Feb 2012 data. The iT s -band photometry presented 
here is based on a weighted mean of the two epochs 1 . 
We note that the changes in JHK values of this order 
correspond to changes of a small fraction of 1 Mj up 
in inferred mass (^0. 1-0.2 Mj up ) regardless of the 
choice of mass-luminosity model, hence there is no 
significant change to the mass estimation presented in 
Kuzuharaet al. (2013). On this topic, it is also worth 
noting that the uncertainties in physical quantities from 
these observational uncertainties are strongly domi- 
nated by model uncertainties, as could also be seen 
in Kuzuhara et al. (2013). For instance, comparing 
the old and the new measurements within the context 
of the Spiegel & Burrows (2012) models, the differ- 
ences in inferred temperature and luminosity are 1% 
and 6%, respectively. By contrast, the difference for 
equal measurements between the Baraffe et al. (2003) 
models and the Spiegel & Burrows (2012) models are 
14% in temperature and 51% in luminosity. This is 
simply because the Spiegel & Burrows (2012) models 
generally predict fainter magnitudes for a given effec- 
tive temperature and luminosity. We use the dominat- 
ing span between the models as the uncertainties when 
we plot the luminosity and temperature in Figs. 2 and 
3. 

Using the full range in luminosity, we show in Fig. 
2 that the mass estimate remains well constrained to 


1 We note here that the newer (May 2013) ft's -band image has in fact 
already been plotted in Kuzuhara et al. (2013). 
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the planetary mass range, irrespective of initial condi- 
tions. The agreement between the photometry of our 
new K' s -band observation with the previous one also 
demonstrates the general robustness of the photome- 
try. However, one issue that should be kept in mind 
is that the 77-band observation displays a puzzling be- 
havior, where the companion appears to have a 2 a dif- 
ferent brightness in the first and second halves of the 
data set, apparently irrespective of calibration issues. 
It may be that something (e.g., temporary cloud cover- 
age) occurred during some part of the observation such 
that it should not be trusted. However, there are no vi- 
sual or quantitative clues in the data to support such 
an assessment. In general, it should be recalled that 
the difference of the two halves of the data set could 
simply be a statistical fluctuation, as 2cr is not a sta- 
tistically unreasonable deviation. We note that, as was 
mentioned by (Kuzuharaet al. 2013), the average fT- 
band photometry of the May 2011 epoch is consistent 
with the photometry from the March 2011 discovery 
epoch, which may support the statistical fluctuation in- 
terpretation. Future observations will be necessary to 
conclusively settle this issue. 

With the CH 4 S and O II 4 L observations taken si- 
multaneously, the methane detection is robust to any 
photometric calibration uncertainties. The absorption 
is deep and appears consistent with late T-dwarfs in the 
field. In order to quantify this, we calculated spectral 
indices from brown dwarfs in the SpeX prism library 2 
using CH 4 S and CH 4 L filter functions. A compari- 
son to GJ 504 b is shown in Fig. 3. It can be seen 
that the methane absorption of the planet is consistent 
with spectral types later than T6. To further check how 
GJ 504 b compares to empirical spectra, we also cal- 
culated JHK photometry from the SpeX library in 
addition to CH 4 S and CH 4 L, and derived \ 2 after 
applying a best-fit scaling factor and penalizing spec- 
tra predicting a CH 4 L detection (Currie et al. 2011). 
Examples from the comparison are shown in Fig. 4. 
This analysis again suggests a late spectral type, with 
T8 providing a better fit than earlier spectral types, al- 
though the best fit still is not impressive, due to the 
excess brightness at /\ -band of GJ 504 b relative to 
field dwarfs already noted in Kuzuhara et al. (2013). 

Methane has otherwise remained elusive in exo- 
planetary atmospheres thus far. While indications of 
methane have been suggested for the transiting planets 


2 Maintained by Adam Burgasser at http://pono.ucsd.edu/ 
~adam/browndwarfs/spexprism 


Table 1: Photometry of GJ 504 b. 


Band 

GJ 504 A 

GJ 504 b 

J 

4. 11 ±0.01 mag 

19.76±0.10 mag 

H 

3.86±0.01 mag 

19.99±0.10 mag a 

K s 

3.81 ±0.01 mag 

19.38±0.11 mag 

CH 4 S 

3.88±0.01 mag 

19.59±0.12 mag 

CH 4 L 

3.85±0.01 mag 

>20.62 mag (3er) 


“There is an apparent systematic difference between the first and second 
halves of the data set, see text for discussion. 



Fig. 2. — Evolution of luminosity with age as a 
function of mass for the coldest (dashed lines) and 
hottest (solid lines) initial conditions of typical warm- 
start models (Spiegel & Burrows 2012). The relatively 
old age of GJ 504 b means that it is only weakly 
sensitive to the initial conditions. Other objects in- 
cluded in the figure are (3 Pic b (Lagrange et al. 2010), 
1RXS J1609 b (Lafreniere et al. 2010), 2M1207 b 
(Chauvin et al. 2005) and the HR 8799 planets 
(Marois et al. 2008, 2010). 
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HD 189733 b and WASP-12 b, these interpretations 
are ambiguous (e.g. Swain et al. 2008; Gibson et al. 
2011; Mandell et al. 2011; Stevenson et al. 2013). 
Also for the case of directly imaged planets, there 
has been a trend of CO/CO2 being the dominant 
carbon-bearing species, rather than CH4 as is expected 
under chemical equilibrium (e.g. Lodders & Fegley 
2002; Konopacky 2013; Lee et al. 2013, and refer- 
ences above). Methane absorption has been identi- 
fied in the atmosphere of HR 8799 b (Barman et al. 
2011), though at modest strength. A stronger case of 
methane absorption has been reported in Ross 458 C 
(Burgasser et al. 2010), which is a very wide (~1000 AU) 
companion to an M-dwarf binary with an inferred mass 
near the deuterium-burning limit. In Kuzuhara et al. 
(2013), we noted that the colors of GJ 504 b are simi- 
lar to those of GJ 758 B (Thalmann et al. 2009), which 
is a similarly cold brown dwarf companion to a Sun- 
like star. Our current detection of methane further 
underlines this similarity, as GJ 758 B also features 
strong methane absorption (Janson et al. 2011). This 
strengthens the indications that GJ 504 b belongs to a 
different class of atmospheres than most of the previ- 
ously imaged planets. The latter share similarities with 
L-type or L/T transition brown dwarfs, while GJ 504 b 
appears to have more in common with the T spectral 
type range. However, like GJ 758 B, it is anoma- 
lously bright in the AT-band range, which could be a 
signature of high metallicity or possibly low surface 
gravity. Indeed, this feature of high A' -band bright- 
ness relative to field brown dwarfs of similar temper- 
ature is a feature that is also present in other objects 
with indications of very low gravity, such as CFBD- 
SIR J2 14947. 2-040308. 9 (Delorme et al. 2012) and 
Ross 458 C (Goldman et al. 2010; Burningham et al. 
2011, e.g.). 

The similarity of GJ 504 b to late T-type ob- 
jects, implying relatively little cloud coverage (al- 
though remnant clouds cannot be excluded, see e.g. 
Burgasser et al. 2010; Morley et al. 2012), is particu- 
larly interesting in the context that thick clouds have 
been inferred even for HR 8799 b, whose temperature 
may be as low as ~800 K (e.g. Barman et al. 2011), 
significantly colder than the temperature of ~ 1 400 K 
at which brown dwarfs lose their thick clouds (e.g. 
Allard et al. 2001). This implies that the transition 
from a more to a less cloudy regime may occur some- 
where in the temperature range of 600-800 K for ob- 
jects in the planetary mass range, and is consistent with 
the predictions of Marley et al. (2012). Future detec- 



T Dwarf Subclass 

Fig. 3. — A comparison of the GJ 504 b methane flux 
ratio limit to the spectral type sequence of field brown 
dwarfs. The green horizontal line marks the measured 
limit of GJ 504 b. The dotted line represents our cal- 
culated indices from the SpeX prism library. 



Fig. 4. — Comparison of GJ 504 b to the calcu- 
lated magnitudes of field brown dwarfs from mesured 
spectra (McLean et al. 2003) in the Spex prism li- 
brary. Only late T-type objects can reproduce the deep 
methane absorption seen in GJ 504 b. 
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tions of planets in this intermediate range would be 
critical to test this hypothesis. The necessary quantity 
of such objects may be delivered by near-future fa- 
cilities such as GPI (Macintosh et al. 2008), SPHERE 
(Beuzitetal. 2008), or CHARIS (Peters et al. 2012). 
Similarly, it will be useful to study planetary mass 
companions to very low-mass stars and brown dwarfs 
(e.g. Chauvin et al. 2005), or even free-floating objects 
in this mass range (e.g. Delorme et al. 2012) in order 
to test for possible similarities or distinctions between 
these different populations, and thus to clarify mutual 
trends with surface gravity, or divergent trends from 
formation scenarios. 
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Fig. 1. — Image showing the methane absorption in GJ 504 b. Left panel: The CH 4 S filter, outside of the methane 
absorption band. Right panel: The CH 4 L filter, inside of the methane absorption band. The planet is bright in CH 4 S 
but completely disappears in CH 4 L, implying strong methane absorption. 
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